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Abstract 

Here we present the two-photon and two-gluon decay widths of the S-wave 
i.VQec,b) and P-wave (xQec,bj) charmonium and bottonium states and the ra- 
diative transition decay widths of cc, bb and cb systems based on Coulomb 
plus power form of the inter-quark potential (CPP U ) with exponent v. The 
Schrodinger equation is solved numerically for different choices of the ex- 
ponent v. We employ the masses of different states and their radial wave 
functions obtained from the study to compute the two-photon and two-gluon 
decay widths and the El and Ml radiative transitions. It is found that the 
quarkonia mass spectra and the El transition can be described by the same 
interquark model potential of the CPP V with v = 1.0 for cc and v = 0.7 
for bb systems, while the Ml transition (at which the spin of the system 
changes) and the decay rates in the annihilation channel of quarkonia are 
better estimated by a shallow potential with v < 1.0. 
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PACS numbers: 12.39.Pn, 12.39.Jh, 14.40.Pq, 13.20.Gd 

1. Introduction 

The renaissance in the hadron spectroscopy particularly in the heavy 
flavour sector associated with the large number of charm and beauty states 
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observed in recent experiments [H, 0, 0, 0] has created renewed interest in 
the study of hadron spectroscopy. Particularly the discovery of i] b (lS) state, 
r] c (2S) state and other orbitally excited states in the charmonia and bottonia 
systems [5[ has stimulated the theoretical phenomenologists to take a new 
look and refine their parameters that describe the nature of quark-antiquark 
interaction at the hadronic scale. The goodness of the spectroscopic parame- 
ters like the interquark potential and its parameters that describe the masses 
of the bound states and the corresponding wave functions obtained from the 
phenomenology, in the descriptions of other properties like the decay (in the 
annihilation channel) and transition properties now become the prospects 
for a detail investigation. With the recent CLEO measurements 0, 5[ of 
the two-photon decay rates of the even-parity, P-wave ++ (xqo) an d 2 ++ 
{XQ2) states and with renewed interest in radiative decays of heavy quarko- 
nium states, it seems appropriate to have another look at the two-photon 
decay of heavy quarkonia from the view point of quark-antiquark interaction 
0, 0j- Though the physics of quarkonium decay seems to be better under- 
stood within the conventional framework of QCD [gj], unlike the two-photon 
width of S-wave (i]q) which can be predicted from the corresponding J/ijj 
and T leptonic widths, the similar prediction for the P-wave (xq) states is 
not conclusive. Similarly the radiative transitions in heavy quarkonia (cc, 



lOj . as they can 



bb and cb) states have drawn much theoretical interest [9|, 
provide direct information on the nature of QQ interaction. 
Most of the existing theoretical values for the decay rates are based on poten- 
tial model calculations that employ different types of interquark potentials 



Hi, \m llS Ud, m Il6|, U3, UM, llS, l20|, l2ll |22|, |23|. We consider the conven- 



tional nonrelativistic formalism for computations of the decay properties of 
the heavy flavour systems. In the traditional non-relativistic bound state 
calculation, the two-photon and two-gluon widths of quarkonium states de 
pend on the £ th derivatives of the radial wave function at the origin 11 



As considered by many authors [24j, |25|, |26|, |27|, |28j], the contribution from 



the radiative corrections to these decays are also incorporated in the present 
study. 

The paper is organized as follows. In Section [2J we describe the phenomeno- 
logical quark-antiquark interaction potential and extract the parameters that 
describe the ground state masses of cc, bb and cb systems. We also compute 
the low lying orbital excited states of these systems. In section [3] we employ 
the spectroscopic parameters of the cc and bb systems to study the two pho- 
ton and two-gluon decay widths. The radiative transitions (El and Ml) of 
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the cc, bb and cb systems are described in section |H In section Owe present, 
analyze and discuss our results to draw important conclusions. 

2. The phenomenology and extraction of the spectroscopic param- 
eters 

For the description of the quarkonium bound states, we adopt the phe 



nomenological Coulomb plus power potential (CPP U ) expressed as [24J, [25 



4 rv 

V(r) = --^-+Ar» (1) 

Here, A is the confinement strength of the potential and a s is the running 
strong coupling constant which is computed as, 

4:71 

a ' ( " 2) = (ii-!„ / )M/.'/A>) (2) 

Where, Uf is the number of flavors, fi is the renormalization scale related to 
the constituent quark mass and A is the QCD scale which is taken as 0.150 
GeV by fixing a s = 0.118 at the Z— boson mass (91 GeV^) [5|. 
The potential parameter, A of EqnJT] is similar to the string strength a of 




Figure 1: Potential strength A (in GeV 1/+1 ) obtained from the ground state spin- average 
mass against the choices of potential index, v (0.1 <v< 2.0) 



the Cornell potential. We particularly chose to vary v in our study as very 
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different interquark potentials can provide fairly good description of the mass 
spectra, while the transitions and other decay properties are very sensitive 
to the choice of interquark potential. Thus the present study on the decay 
properties of heavy flavour mesons based on the CPP V model by varying 
the exponent v (0.1 < v < 2.0) can provide significant understanding of the 
quark-antiquark interaction in the mesonic states while they undergo a tran- 
sition or decay through annihilation channels. The different choices of v here 
then correspond to different potential forms. So, the potential parameter A 
expressed in GeV u+1 can be different for each choices of v. The model poten- 
tial parameter A and the mass parameter of the quark/antiquark (mi, 7712) 
are fixed using the known ground state center of weight (spin average) mass 
and the hyperfine splitting (M-a Sl — Mi So ) of the ground state cc and bb sys- 
tems respectively. The spin average mass for the ground state is computed 
for the different choices of v in the range, 0.1 < v < 2.0. The spin average 
or the center of weight mass, Mew is calculated from the known experimen- 
tal/theoretical values of the pseudoscalar (J = 0) and vector (J = 1) mesonic 
mass as 

£(2J + 1) M nJ 

M n ,cw = J E(2J + 1) (3) 
J 

The Schrodinger equation is numerically solved using the mathematica note- 
book [ver. 3.0] of the Runge-Kutta method [29]. For computing the mass 
difference between different spin degenerate mesonic states, we consider the 
spin dependent part of the usual one gluon exchange potential (OGEP) given 



by |30j, |3l|, |32|, |33|, |34| . Accordingly, the spin- dependent part, VsD(r) contains 
three types of interaction terms, such as the spin-spin, the spin-orbit and the 
tensor part as 



V SD (r) 



V S s(r] 
V T (r) 



3{3 + l) 



S(S + 1) - 



+ V LS (r) L-S) + 



3{S-f){S-f) 



(4) 



The spin-orbit term containing V/^r) and the tensor term containing V^(^) 
describe the fine structure of the meson states, while the spin-spin term 
containing Vss(r) proportional to 2{s q ■ sq) = S(S + 1) — | gives the spin 
singlet-triplet hyperfine splitting. The coefficient of these spin-dependent 
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terms of Eqnf4] can be written in terms of the vector (Vy) and scalar (V$) 
parts of the static potential, V(r) described in Eqn. [T] as [32j 



V LS {r) 



V T (r) = 



Vss(r) 



dV v dVs 



2 m\m<i r V dr 



1 



L/Vv 

6 m!m 2 \ (ir 2 



-VVv 



16 7ra, 



.5(3) ( r -) 



(5) 



(6) 



(7) 



3 m\mi 9 m\m,2 

Here VV is the coulumb part (I s * term) and Vs is the confining part (2 nd 
term) of Eqn. [H The computed masses of the QQi^^Lj) states are listed 
in Table [TJ [2] and [3] for different combinations of Q G b, c. The spectroscopic 
parameters thus correspond to the fitted quark masses, the potential strength 
A, the potential exponent, v and the corresponding radial wave functions. 
The fitted mass parameters are m c = 1.28 GeV/c 2 , = 4.4 GeV/c 2 while 
the potential strength A for each choices of v are shown in Fig. [TJ The £ th 
derivative of the corresponding radial wave functions at the origin obtained 
for each choices of the exponent v are plotted in Fig. [2j 



0.0 0J5 



Figure 2: Square of I th derivative of radial wave function at origin in GeV 2t+3 



3. Two-photon and two-gluon decay widths for quarkonium states 

The extracted model parameters and the radial wavefunctions are be- 
ing employed here to compute the two-photon (r 77 (QQ)) an d two gluon 
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Table 1: Masses (in GeV) of charmonium states (n 2S+1 Lj) in the CPP V mode with 
different choices of v 



v Expt. 



State 




0.1 


0.5 


0.7 


0.9 




1.0 




1.1 




1.5 


2.0 


[5] 


l 3 ^ 


3. 


,076 


3.088 


3.093 


3.097 


3.099 


3. 


,100 


3. 


.106 


3.111 


3.097 


1% 


3. 


,045 


3.008 


2.994 


2.981 


2.976 


2. 


,971 


2. 


.955 


2.940 


2.980 


1 3 P 2 


3. 


,174 


3.357 


3.430 


3.495 


3.524 


3. 


,552 


3. 


.647 


3.737 


3.556 


1 3 P 1 


3. 


,171 


3.347 


3.419 


3.484 


3.514 


3. 


,542 


3. 


.640 


3.736 


3.511 


l 3 Po 


3. 


,167 


3.324 


3.386 


3.441 


3.466 


3. 


,489 


3. 


.570 


3.647 


3.415 


l l Pi 


3. 


,172 


3.350 


3.422 


3.485 


3. 


514 


3, 


,542 


3, 


,636 


3.727 


3.526 


I'D, 


3. 


,215 


3.520 


3.653 


3.774 


3. 


830 


3, 


,830 


4, 


,068 


4.249 




l 3 D 2 


3. 


,214 


3.523 


3.662 


3.792 


3. 


854 


3. 


,914 


4. 


.129 


4.354 




1 3 D 1 


3. 


,214 


3.521 


3.663 


3.796 


3. 


860 


3. 


,923 


4. 


.152 


4.399 


3.770 


l l D 2 


3. 


,214 


3.521 


3.658 


3.785 


3. 


844 


3. 


,901 


4. 


.105 


4.314 




2 3 P 2 


3. 


,121 


3.461 


3.633 


3.803 


3. 


870 


3. 


,970 


4. 


.292 


4.661 


3.929 


2 3 P X 


3. 


,120 


3.455 


3.624 


3.792 


3. 


875 


3. 


,958 


4. 


.277 


4.647 




2 3 P 


3. 


,118 


3.439 


3.599 


3.757 


3. 


835 


3. 


,912 


4. 


.206 


4.541 




2 1 P 1 


3. 


,120 


3.457 


3.626 


3.740 


3. 


877 


3. 


,960 


4. 


.277 


4.643 




2 3 5! 


3. 


,087 


3.326 


3.444 


3.558 


3. 


615 


3. 


,670 


3. 


.881 


4.119 


3.686 


2% 


3. 


,079 


3.290 


3.390 


3.487 


3. 


533 


3, 


,579 


3, 


,749 


3.936 


3.638 
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Table 2: Masses (in GeV) of B c states (n 2S+1 Lj) in the CPP V mode with different choices 
of v 



v Others 



State 


0.1 




0.5 


0.7 


1.0 


1.5 


2.0 


[41 


[42] 


[43] 


[26] 


44] 


1 3 S 1 


6.326 


6 


335 


6.339 


6.343 


6.348 


6.351 


6.416 


6.337 


6.317 


6.332 


6.338 


I 1 S 


6.305 


6 


279 


6.268 


6.256 


6.241 


6.231 


6.380 


6.264 


6.253 


6.270 


6.271 


1 3 P 2 


6.453 


6 


679 


6.771 


6.889 


7.043 


7.156 


6.837 


6.747 


6.743 


6.762 


6.768 


l 3 P l 


6.451 


6 


673 


6.764 


6.882 


7.038 


7.155 


6.772 


6.730 


6.717 


6.734 


6.741 


i 3 Po 


6.448 


6 


656 


6.740 


6.848 


6.989 


7.092 


6.693 


6.700 


6.683 


6.699 


6.706 


i 1 p 1 


6.452 


6 


675 


6.765 


6.882 


7.035 


7.148 


6.775 


6.736 


6.729 


6.749 


6.750 


1 3 D 3 


6.505 


6 


889 


7.059 


7.286 


7.599 


7.840 


7.003 


7.005 


7.007 


7.081 


7.045 


l 3 D 2 


6.504 


6 


891 


7.065 


7.304 


7.643 


7.915 


7.000 


7.012 


7.001 


7.077 


7.041 


1 3 D 1 


6.504 


6 


890 


7.066 


7.308 


7.659 


7.947 


6.959 


7.012 


7.008 


7.720 


7.028 


1 X D* 


6.505 


6 


890 


7.062 


7.297 


7.626 


7.886 


7.001 


7.009 


7.016 


7.079 


7.036 


2 3 P 2 


6.387 


6 


812 


7.027 


7.345 


7.852 


8.313 


7.186 


7.153 


7.134 


7.156 


7.164 


2 3 P l 


6.386 


6 


808 


7.021 


7.337 


7.843 


8.303 


7.136 


7.135 


7.113 


7.126 


7.145 


2 3 P 


6.384 


6 


797 


7.003 


7.309 


7.792 


8.229 


7.081 


7.108 


7.088 


7.091 


7.122 


2 1 P 1 


6.386 


6 


809 


7.022 


7.339 


7.842 


8.301 


7.139 


7.142 


7.124 


7.145 


7.150 




6.343 


6 


640 


6.785 


6.997 


7.326 


7.617 


6.896 


6.899 


6.902 


6.881 


6.887 


2 1 S'o 


6.338 


6 


615 


6.748 


6.939 


7.232 


7.487 


6.875 


6.856 


6.867 


6.835 


6.855 
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Table 3: Masses (in GeV) of bottonia states (n 2S+1 Lj) in the CPP V mode with different 
choices of v 



v Expt. 



State 


0.1 


0.5 


0.7 


0.9 


1.0 


1.5 


2.0 


[5] 


1 S S 1 


9.447 


9.455 


9.458 


9.460 


9.461 


9.466 


9.468 


9.460 


1% 


9.427 


9.404 


9.395 


9.388 


9.385 


9.372 


9.363 


9.389 


1 3 P 2 


9.595 


9.838 


9.938 


10.025 


10.065 


10.231 


10.356 


9.912 


1 3 P, 


9.593 


9.831 


9.929 


10.016 


10.056 


10.223 


10.348 


9.893 


l 3 Po 


9.590 


9.817 


9.909 


9.989 


10.026 


10.179 


10.292 


9.859 




9.594 


9.834 


9.932 


10.018 


10.057 


10.223 


10.346 




1 3 D 3 


9.653 


10.066 


10.250 


10.419 


10.498 


10.842 


11.111 




l 3 D 2 


9.653 


10.067 


10.254 


10.427 


10.509 


10.871 


11.163 


10.162 


1 3 D 1 


9.652 


10.066 


10.253 


10.428 


10.511 


10.881 


11.182 




l x D 2 


9.653 


10.060 


10.252 


10.423 


10.504 


10.859 


11.143 




2 3 P 2 


9.530 


9.984 


10.214 


10.443 


10.556 


11.100 


11.599 


10.269 


2 3 P l 


9.529 


9.970 


10.208 


10.434 


10.547 


11.088 


11.583 


10.255 


2 3 P 


9.528 


9.970 


10.130 


10.413 


10.522 


11.044 


11.518 


10.232 


2 1 P 1 


9.530 


9.981 


10.210 


10.437 


10.549 


11.090 


11.585 




2 3 5i 


9.481 


9.796 


9.950 


10.101 


10.175 


10.524 


10.835 


10.023 


2 1 S'o 


9.476 


9.773 


9.918 


10.057 


10.125 


10.444 


10.724 
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(r gg (QQ)) decay widths. Two photon widths of orbitally excited quarkonium 
states Xqj=o,2 77 are suppressed by the mass of the heavy quark while 
the two-photon decays of the spin one state %Q,i is forbidden by the Landau- 
Yang theorem 35, 36]. Most of the quark model predictions 0, [2l[ 37] for 
the t]q — > 77 width are comparable with the experimental value. However 
the theoretical predictions for the P-wave (xqo,2 -* 77) widths differ largely 
from the experimental observations j^]. This warranted to incorporate con- 
tribution from QCD corrections. Thus with the one-loop QCD radiative 
corrections the decay widths of 1 So (r/g), 3 P Q (xqo) an d 3 P 2 (XQ2) states into 
two photons are computed according to the non relativistic expression given 
by (ill, 0, 0, H M, HQ 



3Qe« 



l^o(0)| 



2 r 



2mf 



1 - 



a s (20 - vr 2 ) 



7T 



(8) 



F~,~,(xqo) 



T 77 (XQ2) 



_ 27Qjal m M XQO \R' 1 (0)\ 2 
2m 5 Q 

4 27Q i e a 2 em M XQ2 \R' 1 (0)\ 



1 + B 



a. 



7T 



15 



2m5 



l + B 2 



7T 



(9) 



(10) 



Here, B = 7r 2 /3 — 28/9 and B 2 = —16/3 are the next to leading order (NLO) 
QCD radiative corrections 40, |45|, |46) and are considered to be the same for 
charmonium and bottonium systems. 

Similarly, the two-gluon decay width of t/q, xqo and \Q2 states are given by 



47 



Tgg(XQ0) = 



a 2 M Vc \R ^ 2 
3rrin 



[1 + C Q (a s /n)] 



ZalM XQQ \R[<^ 



rrif 



[l + C 0Q (a s /n)] 



( 4 \ 3a 2 M Yno \R\(0)\ 2 

r ra (x Q2 ) = l T5 ) s XQ J 11 )l [i + c 2Q (a s M) 



(ii) 

(12) 
(13) 



Here, the quantities in the square brackets are the NLO QCD radiative cor- 
rections |40|, |45|, |46[ and the coefficients, Cq, Cqq and C 2 q for the charmonia 
(Q = c) and bottonia (Q = b) are listed in Table HI 

The value of the radial wave functions |ifeo)| 2 of Eqns. MM are obtained 
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tbp 



Table 4: NLO QCD Correction coefficients for the two-gluon decays 





Q=c 


Q=b 


Cq 


4.8 


4.4 




8.77 


10.0 


C2Q 


-4.827 


-0.1 



from the numerical solution of the schrodinger equation corresponding to the 
CPP V potential model description of the QQ bound states described in sec- 
tion [2j The computed digamma widths of the cc and bb states are listed in 
Table |5] and |6] while the digluon widths of cc and bb states are listed against 
the exponent, v in Table [7] and [8] respectively. 

4. Radiative transitions 

As the mass spectra of the quarkonia states are well known, the investiga- 
tion of radiative transitions become important and such study could help us 
to understand the theory of strong interaction in the nonperturbative regime 
of QCD. The motivation for the present work is to study the radiative transi- 
tions (El and Ml transitions) of the QQ (Q £ b, c) states and to investigate 
the interquark potentials that provides the right description of the quarkonia 
states. The nonrelativistic treatment adopted for the study of these heavy 
flavour mesonic systems allows as to apply the usual multipole expansion in 
electrodynamics to compute the transition between the quarkonia states with 
the emission of a photon. The leading terms in this expansion correspond 
to the El and Ml transitions. The electric dipole term (El) is responsible 
for the transition between the S and P states without changing the spin of 
the quark-antiquark pair, while the magnetic dipole term (Ml) describes the 
transition between S = 1 and S = states without changing relative orbital 
momentum (L) of the quark-antiquark pair. Accordingly, electric transitions 
do not change quark spin. These transitions have AL = ±1 and AS = 0. 
The El radiative transition width between initial state, (n 2s+1 Lj) to final 
state, (n 2s+1 L'j) is given by [481 ]. 



r(< / + 7) = ^2>V + i)sfyM 2 (14) 
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Table 5: Two-photon decay widths (in keV) of cc states with different choices of v 



V 




r 77 (xco) 


r 77 (xc2) 


0.5 


6.375 


2.539 


0.309 


0.7 


7.875 


4.839 


0.529 


1.0 


9.684 


9.357 


1.162 


-i] 




2.36 ±0.35 ±0.22 


0.66 ±0.07 ±0.06 


U 


6.41+JS 


2.397 ±0.4 


0.493 ±0.06 


M 




3.72 


0.49 


[15] 




2.50 


0.28 


\M 




1.39 


0.44 


\XA 




3.34 


0.43 


[isj 




3.78 




W 




1.99 


0.30 


[20] 




1.29 


0.46 


[2JJ 




2.90 


0.50 


[22] 




1.56 


0.56 


[26] 


4.33 






[45J 




3.50 


0.93 


[49] 


4.252 
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Table 6: Two-photon decay widths (in keV) of bb states with different choices of v 



V 


r 77(%) 


r rr (x6o) 


r 7 7(Xb2) 


0.5 


0.410 


0.064 


0.011 


0.7 


0.504 


0.097 


0.016 


1.0 


0.618 


0.185 


0.030 


[13] 




0.080 


0.008 


[15] 


0.460 






[16] 


0.220 ± 0.040 






[12]pert. 


0.30 


0.032 


0.007 


[19]nonpert. 


0.32 


0.094 


0.005 


[20] 


0.214 






[2JJ 


0.350 






[28] 


0.384 ±0.047 






[3JJ 


0.230 






[49] 


0.313 






[50] 


0.466 ±0.101 






[51] 


0.170 






[52] 


0.520 






[53] 


0.560 
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Table 7: Two-gluon decay widths (in MeV) of cc states with different choices of v 



V 




T gg(Xco) 


T gg(Xc2) 


0.5 


32.209 


10.467 


1.169 


0.7 


39.790 


19.949 


2.003 


1.0 


48.927 


38.574 


4.396 


[5] 


26.7 ±3.0 


10.2 ±0.7 


2.03 ±0.12 


[19] pert. 


15.70 


4.68 


1.72 


[19]nonpert. 


10.57 


4.88 


0.69 



Table 8: Two-gluon decay widths (in MeV) of bb states with different choices of v 



V 


T gg(Vb) 


r 99 (xw) 


T gg(Xb2) 


0.5 


11.921 


1.826 


0.285 


0.7 


14.644 


2.745 


0.429 


1.0 


17.945 


5.250 


0.822 


[13] 


12.46 


2.15 


0.22 


[19]pert. 


11.49 


0.96 


0.33 


[19]nonpert. 


12.39 


2.74 


0.25 



Here, a = 1/137, u is the photon energy and (e<g), the mean charge content 
of the QQ system are expressed as 



and 



(eg) 



Mf - M) 
mQeQ ~ m QZQ 



(15) 
(16) 



respectively. Here, Mj and Mf are the initial and final state mass of the 
quarkonia respectively. The statistical factor Sfj = Sf { is given by, 

2 



S 



f = max(£,£)t £ , g £ 



J 1 J 



The overlap integral £j/ is given by 



'-if 



druni(r)u 



nl' (r) 









.T 30 




- Ji It). 



(17) 



(18) 
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Table 9: El transition widths (in keV) of cc states against potential index v 



Initial 


Final 




V 






Others 


State 


State 


0.5 


0.7 




1.0 


[5] 


[55 


1 3 P 2 


l 3 S l 


163.98 


250.22 


383 


48 


406.00 


315 


1 3 P 1 


l 3 S 1 


148.31 


229.62 


361 


21 


320.40 


41 


l 3 Po 


1 3 S 1 


115.53 


173.42 


263 


62 


130.56 


120 


I'Pi 




302.44 


451.50 


671 


05 




482 


1 3 D 3 


1 3 P 2 


137.86 


245.87 


432 


36 




402 


l 3 D 2 


1 3 P 2 


36.23 


68.29 


131 


01 




69.5 


1 3 D 2 


1 3 P! 


127.47 


231.43 


423 


43 




31.3 


l 3 D 1 


1 3 P 2 


3.89 


7.68 


15 


23 




3.88 


I 3 /?! 


1 3 P! 


68.65 


129.97 


245 


70 




99 


l 3 D l 


l 3 Po 


127.71 


240.03 


447 


72 




299 


l l D 2 


1 X P1 


157.16 


285.76 


523 


89 




389 


2 3 P 2 


1 3 S 1 


70.34 


140.10 


252 


68 






2 3 P 2 


2 3 S 1 


36.97 


83.93 


164 


41 






2 3 P X 


l 3 ^ 


66.96 


132.92 


258 


15 






2 3 P l 


2 3 S 1 


32.65 


73.81 


172 


63 






2 3 P 


l 3 S 1 


58.48 


114.28 


216 


53 






2 3 P 


2 3 S 1 


22.63 


49.36 


112 


39 






2 3 P X 


I'So 


122.85 


226.91 


415 


25 






2 x Pi 


2 1 So 


65.04 


146.99 


332 


99 







The El transition rates of charmonia, bottonia and B c systems are listed in 
Table |9j M and EU 

The magnetic dipole transitions (Ml) flip the quark spin, so their ampli- 
tudes are proportional to the quark magnetic moments and thus related in- 
versely to the quark mass. Thus, Ml transitions are weaker than the El tran- 
sitions and this causes difficulties in the experimental observations. Along 
with other exclusive processes, the magnetic dipole (Ml) transitions from the 
spin-triplet S-wave vector (V) state to the spin-singlet S-wave pseudoscalar 
(P) state have also been considered as a valuable testing ground to further 
constrain the phenomenological quark model of hadrons. The Ml transi- 
tions have AL = and the n 2s+1 Lj — > n 2s +1 Ly + 7 transition rate for QQ 
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Table 10: El transition widths (in keV) of bb states against potential index v 



Initial 


Final 




V 




Others 


State 


State 


0.5 


0.7 


1.0 


[5] 


55 


1 3 P 2 


l 3 ^ 


29.06 


44.28 


70.29 


11.88 


31.6 


1 3 P 1 


l 3 5i 


27.58 


42.00 


67.43 


18.91 


27.8 


l 3 Po 


l 3 ^ 


24.75 


37.18 


58.42 


3.24 


22.2 


I'Pi 


1% 


40.29 


60.45 


94.05 






1 3 D 3 


1 3 P 2 


17.28 


41.00 


76.82 




22.6 


l 3 D 2 


1 3 P 2 


4.37 


10.63 


20.60 






1 3 D 2 


1 3 P 


14.32 


34.53 


65.35 






l 3 D l 


1 3 P 2 


0.48 


1.17 


2.32 




0.50 


l 3 Di 


l 3 Pi 


7.86 


19.02 


36.75 




10.7 


1 3 D 1 


l 3 Po 


12.38 


30.04 


58.47 




20.1 


l l D 2 


l^i 


16.84 


44.06 


83.96 






2 3 P 2 




11.05 


20.42 


35.84 




12.7 


2 3 P 2 


2 3 £' 1 


5.54 


13.34 


32.92 




14.5 


2 3 P X 


1 3 S 1 


10.19 


19.92 


34.86 




12.0 


2 3 P X 


2 3 St 


4.43 


12.50 


30.82 




12.4 


2 3 P 


1 3 S 1 


10.19 


14.14 


32.24 




10.9 


2 3 P 


2 3 St 


4.43 


4.44 


25.40 




9.17 


2 x Pi 


1% 


14.39 


25.84 


43.99 
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Table 11: El transition widths (in keV) of B c states against potential index v 



Tnitia 1 

State 


Final 

State 




y 






Others 






7 


1 

X . \J 


44 


26 


41 


42 


l 3 Po 


l 3 ^ 


R? 04 


1 98 1 8 

liU . X (J 


900 07 


83 


107 


109 8 


119 6 

X X Zj . \J 


1 3 P 
X J 1 


13Q 
± <J1 


78 22 


122.73 


1 93 48 


11 


78.9 


67.8 


99 5 

«../ ■ ./ . '.J 


l 3 Po 


l 3 ^ 


67.93 


105.03 


163.13 


60 


67.2 


32.1 


79.2 


I'PI 


1% 


119.47 


184.66 


282.73 






81.8 


56.4 




1 3 P 2 


68.56 


125.28 


225.588 






18.7 


98.7 


1 3 D 2 


1 3 P 2 


17.6 


33.13 


63.37 






4.4 


24.7 


1 3 D 2 


1 3 P 1 


57.06 


105.92 


198.57 






35.8 


88.8 


1 3 D 1 


1 3 P 2 


1.93 


3.71 


7.22 






0.2 


2.7 


1 3 D 1 


l 3 P X 


31.3 


59.38 


113.06 






11.4 


49.3 


l 3 D 1 


l 3 Po 


51.4 


97.17 


183.53 






43.9 


88.6 


l 1 D 2 




73.20 


136.21 


325.46 






120.8 


92.5 


2 3 P 2 


1 3 S 1 


34.78 


68.56 


134.25 






9.1 


25.8 


2 3 P 2 


2 3 St 


17.88 


41.63 


97.24 






91.3 


73.8 


2 3 Pi 


l 3 S 1 


33.89 


66.66 


130.62 






6.3 


22.1 


2 3 Pi 


2 3 St 


16.74 


38.91 


91.62 






53.5 


54.3 


2 3 P 


1 3 S 1 


31.53 


61.18 


118.43 






5.9 


21.9 


2 3 P 


2 3 St 


13.86 


31.36 


73.25 






24.3 


41.2 


2 1 P 1 


l^o 


48.08 


91.81 


174.50 










2 1 P 1 


2% 


24.87 


57.80 


137.46 
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mesonic system is given by 54], |55|, |56 



r(z^/ + 7 ) = ^V^(2J'+l)M 



(19) 



Where 



m Q m c 



(20) 



and, fTi^Q/Q is the mass of the heavy quark/ ant iquark. The statistical factor 



cm _ qM :„ 



J 1 J 



' >v 2 



1 1 

2 2 
S S 



(21) 



The matrix element M contains the integral related to the static term (ij) 
and the term related to the contribution from the two quark confining ex- 
change current (I c ) in the non-relativistic limit as 



Where 



and 



M = h + I c 



/i = </|jo(wr/2)l*> 



m Q 



(f\Ar v \i) 



(22) 



(23) 



(24) 



Here, Ar v is the confining part of the quark- ant iquark potential described by 
CPP V model. 

The magnetic dipole transition rate (allowed Ml transitions) corresponds to 
triplet-singlet transitions between S- wave states, T(n 3 Si — > n 1 ^ + 7) of the 
cc, bb and cb systems are listed in Table [£21 EH and [0] respectively Other 
allowed Ml transitions n 3 Si — > n 1 S$ + 7 with n > n are nonrelativistically 
forbidden 58) and require relativistic treatment which is beyond the scope of 
the present study. The values within the brackets are the results computed 
by incorporating the two-quark exchange current contribution as given by 



56l |57| . Other theoretical model predictions and available experimental 



results are also listed for comparison. 
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Table 12: Ml transition widths (in keV) of cc states against potential index v 



Initial 

State 


Final 






V 




Others 


State 


0.5 


0.7 


0.8 


1.0 


55 26 


57 


\ i S 1 


1% 


1.29(0.47) 


2.41(0.96) 


3.12(1.28) 


4.57(2.01) 


2.0 1.05 


1.05 


2 3 S 1 


2% 


0.12(0.02) 


0.40(0.07) 


0.61(0.10) 


1.39(0.20) 


1.26 


0.14 


3 3 S 1 


3 So 


0.04(0.004) 


0.16(0.009) 


0.28(0.011) 


0.76(0.012) 




0.0124 


Expt 


r(i 3 S! 


-> 1^07) = 


1.21 ±0.37 











Table 13: Ml transition widths (in eV) of bb states against potential index v 



Initial 

State 


Final 

State 






V 






Others 




0.5 


0.7 


0.8 


1.0 


[55] 


[26] 


[57] 


l 3 Si 


l^o 


7.28(5.26) 


13.70(10.11) 


16.47(12.22) 


24.01(18.20) 


8.95 


9.7 


18.9 


2 3 Si 


2% 


0.67(0.41) 


1.80(1.10) 


3.02(1.80) 


6.87(4.01) 


1.51 


1.6 


2.77 


3 3 S 1 




0.19(0.10) 


0.76(0.38) 


1.34(0.64) 


1.98(0.86) 


0.83 


0.9 


1.32 



Table 14: Ml transition widths (in keV) of B c states against potential index v 



Initial 


Final 




V 








Others 




State 


State 


0.5 


0.7 


0.8 


1.0 


26 


[59] [44] 


m 


l 3 ^ 
2 3 S 1 
3 3 S 1 


1% 

2^o 
3 1 Sq 


0.15(0.02) 
0.01(0.01) 
0.005(0.01) 


0.30(0.06) 
0.04(0.06) 
0.02(0.07) 


0.37(0.08) 
0.07(0.15) 
0.03(0.18) 


0.55(0.13) 
0.17(0.48) 
0.09(0.79) 


0.07 
0.03 


0.01 0.08 
0.03 0.01 


0.05 
0.004 
0.08 eV 
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5. Results and Discussions 



The computed masses of charmonium, B c meson and bottonium low lying 
states using the CPP U model are shown in Table HJ [2] and |3] respectively 
with different choices of exponent v. Their respective experimental values 
(PDG average) are also listed for comparison. The best possible choice of 
the exponent iv) for the description of quarkonia (bb & cc) spectra would be 
the one with minimum statistical deviations of the predicted mass spectra 
with the corresponding experimental values. For this we consider about 
ten experimentally known states (n 2S+1 Lj) of cc and bb systems as given in 
Tables [T] and El and compute the root mean square deviations of the predicted 
masses of these states for each choices of v as 



SDmass(u) 



\ 



N ^ 



V r M CPPu {n^Lj) - M exp (n^ Lj ) 



M exp (n 2S +iLj) 



(25) 



The computed values of the SDmass(zz) for cc and bb systems are plotted 
against v in Fig. [3j The figure shows distinct minima around v = 1.0 in 
the case of cc and around v = 0.7 in the case of bb systems. Thus, we 
conclude that charmonia spectra are better described with the cornell-like 
potential with v = 1.0, while the bottonia spectra are better described by a 
relatively flat potential with v = 0.7. In the absence of known experimental 
excited states for B c system, we consider that better predictions for the B c 
spectra must lie within the range of the exponent 0.7 < v < 1.0. And it is 
found that the predicted masses of the B c states for the choices of v in the 
ran ge, 0.7 < v < 0.9 are in good agreement with other theoretical predictions 



M |41|, |42j, |43 



The computed two-photon and two-gluon decay widths of rjg, xqo and XQ2 
states (Q G c, b) using the spectroscopic parameters of CPP U model are com- 
pared with other theoretical predictions in Tables [5j [6l [7] and [8] respectively. 
The experimentally known two-photon and two-gluon widths of the cc states 
are also compared in Table Eland [7] respectively. The two-gluon decay process 
accounts for substantial part of their hadronic decays and hence we have com- 
pared our results of the cc system with their measured hadronic decay widths. 
Though there are many model predictions for the two-photon decay widths of 
Xco and Xc2 states, only very few predictions for the r\ c —> 77 state exist. Our 
predictions for T 77 (cc) with the exponent lying in the range, 0.4 < v < 0.7 
and that for Y gg (cc) in the range, 0.3 < v < 0.5 are in accordance with other 
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Figure 3: rms deviation in mass for quarkonia states 



model predictions. Fig. H] shows the trend lines for the two-photon and two- 
gluon decay widths of the charmonia states with exponent, v. The horizontal 
lines with error bar represents the respective experimental values obtained 
from the PDG data [5|. Similar to the rms deviations computed for the mass 
predictions, the rms deviations of the predicted two-photon widths of the 
charmonia states with their respective experimental values show a minimum 



around v = 0.5 (see Fig. 6(a)) while that in the case of two-gluon widths 



occur around v = 0.4 (See Fig. 6(b)). The shaded regions in Fig. [4] with 



left aligned lines show the neighborhood region of the exponent v at which 
SDmass is minimum for cc states. And the right aligned shaded region show 
the region of v (0.4 < v < 0.7) around which the minimum rms deviations of 
the two-photon widths (SD — 77) occur. The present analysis thus clearly 
indicates that the mass predictions and the annihilation widths studied here 
cannot be explained by any single potential choice. 

In the case of bb systems, Fig. [5] shows the trend lines for the two-photon 
and two-gluon decay widths with v. These decay widths of r/b, \w and Xb2 
states have not been measured. However many other theoretical models pre- 
dicted the two-photon decay of rjb state, while very few model predictions 
exist for Xbo an d Xb2 states. Our two-photon and two-gluon widths of the 
bb systems within the exponent, 0.3 < v < 0.7 are in accordance with other 
theoretical predictions. New experimental results for these decay widths of 
bb system are awaited. Present widths in the range of exponent 0.3 < v < 0.7 



for two-photon and two-gluon decays are marked in Fig. 5(a) as right aligned 
shaded region. 

Further, the present study clearly indicates that the QQ system in its anni- 
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Figure 4: Two-photon and two-gluon decay width of charmonia states against potential 
index v 




Figure 5: Two-photon and two-gluon decay widths of bottonia states against potential 
index v 




Figure 6: rms deviation in two-photon and two-gluon decay of charmonia 
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hilation channel (digamma/digluon decays) interact more weakly compared 
to their bound state spectral excitations. 

The low lying electric dipole transition rates of the cc, bb and cb systems 
are listed in Table [10] and [11] respectively, with other theoretical model 
predictions and with available experimental results. Their trendlines drawn 
against v are shown in Fig. [7] The horizontal lines that cut the El-curves 
correspond to the respective experimental values. The shaded region corre- 
spond to the neighbourhood of v for which the SDmass is minimum. The 
present results for the El-transition rates agree with the experimental data 
[5j in the neighbourhood region of v = 1.0 in the case of charmonium systems 
and in the neighbourhood region of v = 0.7 in the case of bottonium systems. 
These values of v are the same at which SDmass becomes minimum. 
The predicted Ml transition rates with the correction due to the confining 
exchange current contribution for QQ states are drawn in Fig. [8] with re- 
spect to the exponent v. Present result in the case of J/tp — > rjc'-y is in 
agreement with the known experimental result of 1.21 ± 0.37 keV at the 
exponent v w 0.8 with the confining exchange current contribution while it 
agrees with the prediction at the exponent v w 0.5 without the exchange 
current contribution. Our results for cc system in the range, 0.7 < v < 1.0 of 
the exponent, v are in accordance with the values predicted by other models. 
From the present study, the importance of the two -quark confining exchange 
current contribution in the Ml transition widths of cc system becomes very 
clear. In the case of bb and cb states, our predictions with and without the 
exchange current contribution in the range 0.6 < v < 1.0 of the exponent, v 
are in accordance with other model predictions. However lack of experimen- 
tal data for these decay widths does not allow to draw conclusion in favour 
of a particular choice of the exponent v. 

To summarize, we find that the description of the quarkonia mass spectra 
and the El transition can be described by the same interquark model poten- 
tial of the CPP V with v = 1.0 for cc and v = 0.7 for bb systems, while the 
Ml transition (at which the spin of the system changes) and the decay rates 
in the annihilation channel of quarkonia are better estimated by a shallow 
potential with v < 1.0. We look forward to future experimental data related 
to the decay properties of cb and bb systems. 
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Figure 7: El transition decay widths for quarkonium states against potential index v 




Figure 8: Ml transition decay widths for quarkonium states against potential index v 
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